Complete mitochondrial genome of Thuja sutchuenensis and its implications on evolutionary analysis of complex mitogenome architecture in Cupressaceae.

BACKGROUND: The complex physical structure and abundant repeat sequences make it difficult to assemble the mitogenomes of seed plants, especially gymnosperms. Only approximately 33 mitogenomes of gymnosperms have been reported. However, as the most widely distributed and the second largest family among gymnosperms, Cupressaceae has only six assembled mitogenomes, including five draft mitogenomes and one complete mitogenome, which has greatly hindered the understanding of mitogenome evolution within this large family, even gymnosperms. RESULTS: In this study, we assembled and validated the complete mitogenome of Thuja sutchuenensis, with a size of 2.4 Mb. Multiple sequence units constituted its complex structure, which can be reduced to three linear contigs and one small circular contig. The analysis of repeat sequences indicated that the numbers of simple sequence repeats increased during the evolutionary history of gymnosperms, and the mitogenome of Thuja sutchuenensis harboured abundant extra-long repeats (more than 5 kb). Additionally, the longest repeat sequence identified in these seven gymnosperms also came from the mitogenome of Thuja sutchuenensis, with a length of up to 47 kb. The analysis of colinear blocks and gene clusters both revealed that the orders of mitochondrial genes within gymnosperms was not conserved. The comparative analysis showed that only four tRNAs were shared by seven gymnosperms, namely, trnD-GUC, trnE-UUC, trnI-CAU and trnY-GUA. Furthermore, four genes have undergone potential positive selection in most gymnosperm species, namely, atp8, ccmB, mttB and sdh4. CONCLUSION: We successfully assembled the second complete mitogenome within Cupressaceae and verified that it consisted of multiple sequence units. Our study also indicated that abundant long repeats may contribute to the generation of the complex conformation of the mitogenome of Thuja sutchuenensis. The investigation of Thuja sutchuenensis's mitogenome in our study provides new insight into further understanding the complex mitogenome architecture within gymnosperms.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization.

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

Comparative analyses of six complete chloroplast genomes from the genus Cupressus and Juniperus (Cupressaceae).

Cupressaceae is a conifer family distributed around the world. Cupressus and Juniperus are the main genera of the Cupressaceae family and have important medicinal value. This leads to confusion between Cupressus and Juniperus due to similar morphologies. Here, the complete cp genomes of two Cupressus (C. duclouxiana and C. funebri) and four Juniperus (J. chinensis, J. gaussenii J. pingii and J. procumbens) were sequenced. The results revealed that the length of the cp genomes ranged from 126,996 bp to 129,959 bp, with 119 genes comprising 82 protein-coding genes, 33 transfer RNAs and 4 ribosomal RNAs. All chloroplast genomes of Cupressus and Juniperus lost whole IR regions, which is consistent with gymnosperm cp genome studies. In addition, the number of SSRs per species ranged from 54 to 73 and was dominated by mononucleotide repeats. In the six cp genomes of Cupressus and Juniperus, five highly divergent regions, including accD, accD-rpl2, ycf1, ycf2 and rrn23-rrn4.5, can be used as DNA barcodes of interspecific relationships and potential genetic markers. We compared the gene selection pressures (C. chengiana as reference species), and 6 genes underwent positive selection, the majority of which were related to photosynthesis. Phylogenetic results showed that the monophyly of Cupressus and Juniperus supported most bootstrap support. Cupressus funebris and J. chinensis were resolved to be early diverging species within Cupressus and Juniperus, and the two genera were sister groups to each other. This research revealed a new understanding of the structural pluralism and phylogenetic relationships of Cupressaceae cp genomes. These results will facilitate comprehension of the complexity and diversity of conifer cp genomes.

Building a reference transcriptome for Juniperus squamata (Cupressaceae) based on single-molecule real-time sequencing.

OBJECTIVES: Cupressaceae is the second largest family of coniferous trees (Coniferopsida) with important economic and ecological values. However, like other conifers, the members of Cupressaceae have extremely large genome (> 8 gigabytes), which limited the researches of these taxa. A high-quality transcriptome is an important resource for gene discovery and annotation for non-model organisms. DATA DESCRIPTION: Juniperus squamata, a tetraploid species which is widely distributed in Asian mountains, represents the largest genus, Juniperus, in Cupressaceae. Single-molecule real-time sequencing was used to obtain full-length transcriptome of Juniperus squamata. The full-length transcriptome was corrected with Illumina RNA-seq data from the same individual. A total of 47,860 non-redundant full-length transcripts, N50 of which was 2839, were obtained. A total of 57,393 simple sequence repeats were identified and 268,854 open reading frames were predicted for Juniperus squamata. A BLAST alignment against non-redundant protein database was conducted and 10,818 sequences were annotated in Gene Ontology database. InterPro analysis shows that 30,403 sequences have been functionally characterized against its member database. This data presents the first comprehensive transcriptome characterization of Juniperus species, and provides an important reference for researches on the genomics and evolutionary history of Cupressaceae plants and conifers in the future.

Selection and Verification of Appropriate Reference Genes for Expression Normalization in Cryptomeria fortunei under Abiotic Stress and Hormone Treatments.

Cryptomeria fortunei has become one of the main timber afforestation species in subtropical high-altitude areas of China due to its fast growth, good material quality, and strong adaptability, showing broad application prospects. Quantitative real-time PCR (qRT-PCR) is the most accurate and widely used gene expression evaluation technique, and selecting appropriate reference genes (RGs) is essential for normalizing qRT-PCR results. However, suitable RGs for gene expression normalization in C. fortunei have not been reported. Here, we tested the expression stability for 12 RGs in C. fortunei under various experimental conditions (simulated abiotic stresses (cold, heat, drought, and salinity) and hormone treatments (methyl jasmonate, abscisic acid, salicylic acid, and gibberellin) and in different tissues (stems, tender needles, needles, cones, and seeds) using four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper). Then, geometric mean rankings from these algorithms and the RefFinder program were used to comprehensively evaluate RG stability. The results indicated CYP, actin, UBC, and 18S as good choices for studying C. fortunei gene expression. qRT-PCR analysis of the expression patterns of three target genes (CAT and MAPK1/6) further verified that the selected RGs were suitable for gene expression normalization. This study provides an important basis for C. fortunei gene expression standardization and quantification.

Development of polymorphic microsatellite markers for the Tertiary relict tree species Taiwania cryptomerioides (Cupressaceae) in East Asia.

Taiwania cryptomerioides Hayata is an endangered relict tree species which is endemic to mainland China, Taiwan, Myanmar, and northern Vietnam. It is an economically important tree species and has been used for reforestation in mountain areas of mainland China and Taiwan. In order to investigate its genetic diversity for conservation and restoration, we developed and characterized 15 nuclear microsatellite markers based on next-generation sequencing data. A total of 100 microsatellite primer pairs were initially designed and tested based on the restriction-site associated DNA sequencing data. 60 of 100 loci (60%) were successfully amplified, of which 42 loci exhibited polymorphism. Fifteen polymorphic microsatellite loci with clear peaks were selected for further analyses in four T. cryptomerioides populations sampled from China (Hubei, Fujian, Guizhou, and Yunnan). The number of alleles per locus ranged from 2 to 24, and the levels of observed and expected heterozygosity ranged from 0.000 to 0.950 and from 0.000 to 0.860, respectively. This set of microsatellite markers will be useful for future population genetic studies of T. cryptomerioides in East Asia.

Novel 28 microsatellite loci using high-throughput sequencing for an endangered species on Metasequoia glyptostroboides (Cupressaceae).

Metasequoia glyptostroboides is a living fossil and an endangered species listed in the International Union for Conservation of Nature (IUCN). Distinguishing the genotypes of all wild individuals of M. glyptostroboides is important to delimit management units and key germplasm resources. We characterized 28 novel polymorphic microsatellite loci using a streptavidin-biotin microsatellite-enriched library and Illumina high-throughput sequencing. Characteristics of each locus were tested using 140 individuals collected from five natural populations of M. glyptostroboides. The number of alleles per locus ranged from 3 to 20, with a mean number of about 8 alleles. The observed and expected heterozygosities in each population ranged from 0.0000 to 1.0000 and from 0.0000 to 0.8958, respectively. Four to nine loci were cross-amplified successfully in seven species of Cupressaceae. The novel SSR markers will provide a toolkit for DNA identification of all of the extant wild individuals guiding further conservation efforts of M. glyptostroboides.

Differential Gene Profiling of the Heartwood Formation Process in Taiwania cryptomerioides Hayata Xylem Tissues.

Taiwania (Taiwania cryptomerioides) is an important tree species in Taiwan because of the excellent properties of its wood and fascinating color qualities of its heartwood (HW), as well as the bioactive compounds therein. However, limited information is available as to the HW formation of this species. The objective of this research is to analyze the differentially expressed genes (DEGs) during the HW formation process from specific Taiwania xylem tissues, and to obtain genes that might be closely associated with this process. The results indicated that our analyses have captured DEGs representative to the HW formation process of Taiwania. DEGs related to the terpenoid biosynthesis pathway were all up-regulated in the transition zone (TZ) to support the biosynthesis and accumulation of terpenoids. Many DEGs related to lignin biosynthesis, and two DEGs related to pinoresinol reductase (PrR)/pinoresinol lariciresinol reductase (PLR), were up-regulated in TZ. These DEGs together are likely involved in providing the precursors for the subsequent lignan biosynthesis. Several transcription factor-, nuclease-, and protease-encoding DEGs were also highly expressed in TZ, and these DEGs might be involved in the regulation of secondary metabolite biosynthesis and the autolysis of the cellular components of ray parenchyma cells in TZ. These results provide further insights into the process of HW formation in Taiwania.

Phylogenetically distant group of terpene synthases participates in cadinene and cedrane-type sesquiterpenes accumulation in Taiwania cryptomerioides.

Along with the species evolution, plants have evolved ways to produce a different collection of terpenoids to accommodate its biotic and abiotic environment, and terpene synthase (TPS) is one of the major contributors to various terpene compounds. The timber of a monotypic and relictual conifer species of Cupressace, Taiwania cryptomerioides, has excellent durability, and one of the essential factors for Taiwania to resist decay and insect pests is sesquiterpene. Compared to other conifers, Taiwania has much higher abundance of cadinene-type sesquiterpenes, and the presence of cedrene-type sesquiterpenes. To understand sesquiterpene biosynthesis in Taiwania, we functionally characterized 10 T. cryptomerioides TPSs (TcTPSs) in vivo or in planta, which could catalyze sesquiterpene formation and potentially are involved in biosynthesis of diverse sesquiterpenoids in Taiwania. The distant phylogenetic relationship and the intron loss event of TcTPSs correlate to the differentiation of chemical profile Taiwania compared to other conifers. Furthermore, we identified TcTPS3 and TcTPS12 as δ-cadinene synthase, and TcTPS6 as cedrol synthase, which demonstrates the important contributions of dynamic evolution in TPSs to the chemical diversity in plants. Combining with functional characterization and comparison of catalytic residues, we conclude at least three catalytic routes for sesquiterpene biosynthesis in this species, and the skeleton diversity has been expended in T. cryptomeriodes.

Genome-Wide Variant Identification and High-Density Genetic Map Construction Using RADseq for Platycladus orientalis (Cupressaceae).

Platycladus orientalis is an ecologically important native conifer in Northern China and exotic species in many parts of the world; however, knowledge about the species' genetics and genome are very limited. The availability of well-developed battery of genetic markers, with large genome coverage, is a prerequisite for the species genetic dissection of adaptive attributes and efficient selective breeding. Here, we present a genome-wide genotyping method with double-digestion restriction site associated DNA sequencing (ddRAD-seq) that is effective in generating large number of Mendelian markers for genome mapping and other genetic applications. Using 139 megagametophytes collected from a single mother tree, we assembled 397,226 loci, of which 108,683 (27.4%) were polymorphic. After stringent filtering for 1:1 segregation ratio and missing rate of <20%, the remaining 23,926 loci (22% of the polymorphic loci) were ordered into 11 linkage groups (LGs) and distributed across 7,559 unique positions, with a total map length of 1,443 cM and an average spacing of 0.2 cM between adjacent unique positions. The 11 LGs correspond to the species' 11 haploid genome chromosome number. This genetic map is among few high-density maps available for conifers to date, and represents the first genetic map for P. orientalis The information generated serves as a solid foundation not only for marker-assisted breeding efforts, but also for comparative conifer genomic studies.

Biochemical characterization of diterpene synthases of Taiwania cryptomerioides expands the known functional space of specialized diterpene metabolism in gymnosperms.

Taiwania cryptomerioides is a monotypic gymnosperm species, valued for the high decay resistance of its wood. This durability has been attributed to the abundance of terpenoids, especially the major diterpenoid metabolite ferruginol, with antifungal and antitermite activity. Specialized diterpenoid metabolism in gymnosperms primarily recruits bifunctional class-I/II diterpene synthases (diTPSs), whereas monofunctional class-II and class-I enzymes operate in angiosperms. In this study, we identified a previously unrecognized group of monofunctional diTPSs in T. cryptomerioides, which suggests a distinct evolutionary divergence of the diTPS family in this species. Specifically, five monofunctional diTPS functions not previously observed in gymnosperms were characterized, including monofunctional class-II enzymes forming labda-13-en-8-ol diphosphate (LPP, TcCPS2) and (+)-copalyl diphosphate (CPP, TcCPS4), and three class-I diTPSs producing biformene (TcKSL1), levopimaradiene (TcKSL3) and phyllocladanol (TcKSL5), respectively. Methyl jasmonate (MeJA) elicited the accumulation of levopimaradiene and the corresponding biosynthetic diTPS genes, TcCPS4 and TcKSL3, is consistent with a possible role in plant defense. Furthermore, TcCPS4 and TcKSL3 are likely to contribute to abietatriene biosynthesis via levopimaradiene as an intermediate in ferruginol biosynthesis in Taiwania. In conclusion, this study provides deeper insight into the functional landscape and molecular evolution of specialized diterpenoid metabolism in gymnosperms as a basis to better understand the role of these metabolites in tree chemical defense.

De novo assembly of the Platycladus orientalis (L.) Franco transcriptome provides insight into the development and pollination mechanism of female cone based on RNA-Seq data.

For seed-bearing plants, the basis of seed and fruit formation is pollination. The normal progression of pollination is through advances in continuous signal exchange and material transfer, which occur mainly in female reproductive organs; thus, the molecular mechanism of development in female reproductive organs is vital for understanding the principle of pollination. However, molecular biology studies on the development of female cones related to pollination are rare and unclear in gymnosperms, especially in Cupressaceae. In this study, Platycladus orientalis, a monotypic genus within Cupressaceae, was chosen to examine female cone transcriptomes at pre-pollination and pollination stages by Illumina paired-end sequencing technology to de novo sequence six libraries with 3 biological replicates. These libraries were used to construct a P. orientalis transcriptome database containing 71,669 unigenes (4,963 upregulated unigenes and 11,747 downregulated unigenes at the pollination stage) for subsequent analysis. Based on the annotations and expression levels, the functions of differentially expressed unigenes and enriched pathways between the developmental processes of female cones were analysed to detail the preliminary development and pollination mechanism of the female cone. Targeted investigations were specifically performed to determine the elementary mechanism of secretion and functioning of the pollination drop, a vital ovule secretion at the pollination stage. Ultimately, the expression of 15 unigenes selected between two stages were further assessed and confirmed using qRT-PCR, which demonstrated reliable data and significant differences in the expression profiles of key genes. As one of the largest available transcriptomic resources of this species, the database is constructed to prospectively adapt to the physiological and genomic data of woody plants. This work provided the first transcriptome profile of P. orientalis female cones at different developmental stages, and will promote the illumination of the pollination mechanism of P. orientalis, and will serve as the basis for in-depth genomic study in the Cupressaceae family. This initiative will arouse the interest and attention of scholars and pave the way for future studies.

Genetic diversity and population structure of two endemic Cupressus (Cupressaceae) species on the Qinghai-Tibetan plateau.

Cupressus gigantea and C. torulosa are ecologically and economically important endemic species of the conifer family Cupressaceae on the Qinghai-Tibetan plateau. C. gigantea was previously classified as a subspecies of C. torulosa because of their similar morphological characteristics and close distribution. In this study, 401 individuals were sampled from 16 populations of the two Cupressus species. The specimens were genotyped using 10 polymorphic microsatellite loci through fluorescence polymerase chain reaction (PCR). The genetic diversity of C. gigantea and C. torulosa populations was generally low, with the highest genetic diversity detected in the population LLS of C. gigantea. Distance-based phylogenetic and principal co-ordinates analyses indicated a clear genetic structures for the 16 populations of the two Cupressus species. Moreover, Mantel test results showed indistinctive correlations between population-pairwise Fst values and geographic distances, as well as between genetic distances and geographic distances in C. gigantea and C. torulosa, respectively. AMOVA suggested that genetic variation mostly resided within populations. Sixteen naturalpopulations were evidently clustered into two major groups in the constructed neighbour-joining tree. The results demonstrated that C. gigantea and C. torulosa are different Cupressus species. The genetic information provided important theoretical references for conservation and management of the two endangered Cupressus species.

A transcriptome-based resolution for a key taxonomic controversy in Cupressaceae.

Background and Aims: Rapid evolutionary divergence and reticulate evolution may result in phylogenetic relationships that are difficult to resolve using small nucleotide sequence data sets. Next-generation sequencing methods can generate larger data sets that are better suited to solving these puzzles. One major and long-standing controversy in conifers concerns generic relationships within the subfamily Cupressoideae (105 species, approx. 1/6 of all conifers) of Cupressaceae, in particular the relationship between Juniperus, Cupressus and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) clade. Here we attempt to resolve this question using transcriptome-derived data. Methods: Transcriptome sequences of 20 species from Cupressoideae were collected. Using MarkerMiner, single-copy nuclear (SCN) genes were extracted. These were applied to estimate phylogenies based on concatenated data, species trees and a phylogenetic network. We further examined the effect of alternative backbone topologies on downstream analyses, including biogeographic inference and dating analysis. Results: Based on the 73 SCN genes (>200 000 bp total alignment length) we considered, all tree-building methods lent strong support for the relationship (HCX, (Juniperus, Cupressus)); however, strongly supported conflicts among individual gene trees were also detected. Molecular dating suggests that these three lineages shared a most recent common ancestor approx. 60 million years ago (Mya), and that Juniperus and Cupressus diverged about 56 Mya. Ancestral area reconstructions (AARs) suggest an Asian origin for the entire clade, with subsequent dispersal to North America, Europe and Africa. Conclusions: Our analysis of SCN genes resolves a controversial phylogenetic relationship in the Cupressoideae, a major clade of conifers, and suggests that rapid evolutionary divergence and incomplete lineage sorting probably acted together as the source for conflicting phylogenetic inferences between gene trees and between our robust results and recently published studies. Our updated backbone topology has not substantially altered molecular dating estimates relative to previous studies; however, application of the latest AAR approaches has yielded a clearer picture of the biogeographic history of Cupressoideae.

Comparison of Reliable Reference Genes Following Different Hormone Treatments by Various Algorithms for qRT-PCR Analysis of Metasequoia.

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is the most sensitive technique for evaluating gene expression levels. Choosing appropriate reference genes for normalizing target gene expression is important for verifying expression changes. Metasequoia is a high-quality and economically important wood species. However, few systematic studies have examined reference genes in Metasequoia. Here, the expression stability of 14 candidate reference genes in different tissues and following different hormone treatments were analyzed using six algorithms. Candidate reference genes were used to normalize the expression pattern of FLOWERING LOCUS T and pyrabactin resistance-like 8. Analysis using the GrayNorm algorithm showed that ACT2 (Actin 2), HIS (histone superfamily protein H3) and TATA (TATA binding protein) were stably expressed in different tissues. ACT2, EF1α (elongation factor-1 alpha) and HIS were optimal for leaves treated with the flowering induction hormone solution, while Cpn60ß (60-kDa chaperonin ß-subunit), GAPDH (glyceraldehyde-3-phosphate dehydrogenase) and HIS were the best reference genes for treated buds. EF1α, HIS and TATA were useful reference genes for accurate normalization in abscisic acid-response signaling. Our results emphasize the importance of validating reference genes for qRT-PCR analysis in Metasequoia. To avoid errors, suitable reference genes should be used for different tissues and hormone treatments to increase normalization accuracy. Our study provides a foundation for reference gene normalization when analyzing gene expression in Metasequoia.

Phylogenomic evidence for ancient recombination between plastid genomes of the Cupressus-Juniperus-Xanthocyparis complex (Cupressaceae).

BACKGROUND: Phylogenetic relationships among Eastern Hemisphere cypresses, Western Hemisphere cypresses, junipers, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 11 new plastid genomes (plastomes) from the five variously described genera in this complex (Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis) and compared them with additional plastomes from diverse members of Cupressaceae. RESULTS: Phylogenetic analysis of protein-coding genes recovered a topology in which Juniperus is sister to Cupressus, whereas a tree based on whole plastomes indicated that the Callitropsis-Hesperocyparis-Xanthocyparis (CaHX) clade is sister to Cupressus. A sliding window analysis of site-specific phylogenetic support identified a ~ 15 kb region, spanning the genes ycf1 and ycf2, which harbored an anomalous signal relative to the rest of the genome. After excluding these genes, trees based on the remainder of the genes and genome consistently recovered a topology grouping the CaHX clade and Cupressus with strong bootstrap support. In contrast, trees based on the ycf1 and ycf2 region strongly supported a sister relationship between Cupressus and Juniperus. CONCLUSIONS: These results demonstrate that standard phylogenomic analyses can result in strongly supported but conflicting trees. We suggest that the conflicting plastomic signals result from an ancient introgression event involving ycf1 and ycf2 that occurred in an ancestor of this species complex. The introgression event was facilitated by plastomic recombination in an ancestral heteroplasmic individual carrying distinct plastid haplotypes, offering further evidence that recombination occurs between plastomes. Finally, we provide strong support for previous proposals to recognize five genera in this species complex: Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis.

Transcriptome differences between 20- and 3, 000-year-old Platycladus orientalis reveal that ROS are involved in senescence regulation

Background: Platycladus orientalis has an extremely long life span of several thousands of years, attracting great interests in the mechanisms involved in such successful senescence regulation and resistance at physiological and molecular levels. Results: The levels of reactive oxygen species (ROS) were higher in 3,000-year-old than in 20-year-old P. orientalis, and the activities of GR and GSH demonstrated the same trend. We produced and analyzed massive sequence information from pooled samples of P. orientalis through transcriptome sequencing, which generated 51,664 unigenes with an average length of 475 bp. We then used RNA-seq analysis to obtain a high-resolution age­course profile of gene expression in 20- and 3,000-year-old P. orientalis individuals. Totally, 106 differentially expressed genes were obtained, of which 47 genes were downregulated and 59 upregulated in the old tree. These genes were involved in transcription factors, hormone-related responses, ROS scavengers, senescence-related responses, stress response, and defense and possibly play crucial roles in tackling various stresses in the 3,000-year-old P. orientalis during its life time. The expression patterns of genes related to ROS homeostasis further indicated that the high ability of ROS scavenging could be helpful for the 3,000-year-old P. orientalis to resist senescence. Conclusions: This study provides a foundation for the elucidation of senescence resistance through molecular studies and the discovery of useful genes in P. orientalis.

Aridity drove the evolution of extreme embolism resistance and the radiation of conifer genus Callitris.

Xylem vulnerability to embolism is emerging as a major factor in drought-induced tree mortality events across the globe. However, we lack understanding of how and to what extent climate has shaped vascular properties or functions. We investigated the evolution of xylem hydraulic function and diversification patterns in Australia's most successful gymnosperm clade, Callitris, the world's most drought-resistant conifers. For all 23 species in this group, we measured embolism resistance (P50 ), xylem specific hydraulic conductivity (Ks ), wood density, and tracheary element size from natural populations. We investigated whether hydraulic traits variation linked with climate and the diversification of this clade using a time-calibrated phylogeny. Embolism resistance varied widely across the Callitris clade (P50 : -3.8 to -18.8 MPa), and was significantly related to water scarcity, as was tracheid diameter. We found no evidence of a safety-efficiency tradeoff; Ks and wood density were not related to rainfall. Callitris diversification coincides with the onset of aridity in Australia since the early Oligocene. Our results highlight the evolutionary lability of xylem traits with climate, and the leading role of aridity in the diversification of conifers. The uncoupling of safety from other xylem functions allowed Callitris to evolve extreme embolism resistance and diversify into xeric environments.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae).

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

Genetic evaluation of the breeding population of a valuable reforestation conifer Platycladus orientalis (Cupressaceae).

Platycladus orientalis, a widespread conifer with long lifespan and significant adaptability. It is much used in reforestation in north China and commonly planted in central Asia. With the increasing demand for plantation forest in central to north China, breeding programs are progressively established for this species. Efficient use of breeding resources requires good understanding of the genetic value of the founder breeding materials. This study investigated the distribution of genetic variation in 192 elite trees collected for the breeding program for the central range of the species. We developed first set of 27 polymorphic EST-derived SSR loci for the species from transcriptome/genome data. After examination of amplification quality, 10 loci were used to evaluate the genetic variation in the breeding population. We found moderate genetic diversity (average He = 0.348) and low population differentiation (Fst = 0.011). Extensive admixture and no significant geographic population structure characterized this set of collections. Our analyses of the diversity and population structure are important steps toward a long-term sustainable deployment of the species and provide valuable genetic information for conservation and breeding applications.